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ABSTRACT 

We study the effect of inhomogeneities in the matter distribution of the universe 
on the Faraday rotation of light from distant QSOs and derive new limits on the 
cosmological magnetic field. The matter distribution in the Universe is far from being 
homogeneous and, for the redshifts of interest to rotation measures (RM), it is well 
described by the observed Ly-a forest. We use a log- normal distribution to model the 
Ly-a forest and assume that a cosmological magnetic field is frozen into the plasma 
and is therefore a function of the density inhomogeneities. The Ly-a forest results are 
much less sensitive to the cosmological magnetic field coherence length than those for a 
homogeneous universe and show an increase in the magnitude of the expected RM for 
a given field by over an order of magnitude. The forest also introduces a large scatter 
in RM for different lines-of-sight with a highly non-gaussian tail that renders the 
variance and the mean RM impractical for setting limits. The median\KM.\ is a better 
statistical indicator which we use to derive the following limits using the observed RM 
for QSOs between z = and z = 2.5. We set Qbh 2 = 0.02 and get for cosmological 
fields coherent accross the present horizon, B H -i <^ 10~ 9 G in the case of a Ly-a forest 
which is stronger than the limit for a homogeneous universe, -B^-i ^ 2 x 10~ 8 G; 

while for 50 Mpc coherence length, the inhomogeneous case gives -BsoMpc ^ 6 x 10 -9 
G while the homogeneous limit is B^ 0M <^ 10~ 7 G and for coherence length equal to 
the Jeans length, B\ r £ 10~ 8 G for the Ly-a case while B% £ 10 -6 G. 



1. Introduction 

The nature and structure of extragalactic magnetic fields are currently unknown. Magnetic 
fields in clusters of galaxies and in bridges between clusters have been observed (Kim et al. 1989; 
Vallee 1990; Kronberg 1994 and references therein) but fields in cosmological scales are yet to be 
established. Cosmological fields produced before galaxy formation would be present throughtout 
the Universe today and could be detected by Faraday rotation observations. A number of models 
for the generation of such primordial fields have been suggested (e.g., Harrison 1969, 1970, 1973; 
Hogan 1983; Turner & Widrow 1988; Vachaspati 1991; Quashnock, Loeb, & Spergel 1989; Ratra 
1992; Enqvist & Olesen 1993; Dolgov & Silk 1993; Cheng & Olinto 1994; Kibble & Vilenkin 1995; 
Gasperini, Giovannini, & Veneziano 199, Lemoine & Lemoine 1995; Baym, Bodecker, & McLerran 
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1996; Sigl, Olinto, & Jedamzik 1997; Joyce & Shaposhnikov 1997; Giovannini & Shaposhnikov 
1997) with hopes of seeding magnetic fields in galaxies. The evolution of these fields has been 
studied before (Jedamzik, Katalinic, & Olinto 1998; Brandenburg, Enqvist, & Olesen 1997) and 
around recombination (Wasserman 1978; Kim, Olinto, h Rosner 1996) but theoretical expectations 
for the present structure of these fields has not been contrasted with observations. Here we take 
a first step in this direction by considering the effect of densities inhomogeneities in the present 
Universe on Faraday rotation measures and find much stronger limits on the primordial field than 
the previously assumed case of a homogeneous density distribution. 

As early researchers realized (Rees & Reinhardt 1972; Nelson 1973; Kronberg h Simard- 
Normandin 1976), the most effective method for studying large scale magnetic fields is the 
detection of Faraday rotation from extragalactic sources. Faraday rotation occurs when polarized 
electromagnetic radiation travels through a magnetized medium. The polarization vector rotates 
an angle, ip, proportional to the magnetic field along the line-of-sight, B\\(l), and varies with the 
observed wavelength as \~ bs . The rotation measure is defined as, 



where n e is the number density of electrons (with charge e and mass m e ), A(Z) is the wavelength at 
a physical position / along the line-of-sight to a source at physical distance l s , and c is the speed 



This method has been quite sucessful in determining the structure of the Galactic magnetic 
field by observing both RM and the dispersion measure of pulsar radio emission. With the 
dispersion measure data, the electron number density can be estimated and a determined. For 
extragalactic sources, no dispersion measure data is attainable and the electron density has to be 
modeled theoretically. 

Prior studies of extragalactic RM have generally used oversimplified models for the electron 
density distribution given what we presently know about the matter distribution in the Universe. 
In fact, the electron density is usually assumed to be a homogeneous distribution that only scales 
with redshift, z, as n e oc (1 + z) 3 . Instead, it has been clear for a number of years that the matter 
in the intergalactic medium is not homogeneously distributed. A new picture of the intergalactic 
medium has emerged in which the baryon fluctuations trace the dark matter fluctuations on 
scales larger than the Jeans length (Cen et al. 1994). This picture fits well with high resolution 
spectroscopic observations of the ubiquitous high-redshift neutral hydrogen absorption clouds 
known as the Lyman-a forest (Kirkman & Tylter 1997, Hu et al. 1995, Lu et al. 1996) and with 
N-body+hydro simulations of large-scale structure formation (Rauch et al. 1997, Zhang et al. 



In this Letter, we study the RM of distant QSOs assuming that the electron distribution 
follows the observed Lyman-a forest distribution. We show that not only the variance and 
mean RM for different lines-of-sight can differ by more than an order of magnitude from the 
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case of a homogeneous universe but their scatter is too large to determine a significant limit. 
The dependence on the cosmological magnetic field coherence scale is almost washed out 
becoming sub dominant to the density fluctuations scale. To determine a limit, we use instead the 
median\RM\ and find that for Qbh 2 = 0.02 and Qm = 1 the cosmological field is B <^ 10~ 9 — 10~ 8 
G for coherence lengths varying from the Jeans length to the horizon scale, Hq 1 . 

In fact, for cosmological fields coherent accross the present horizon, the case of a Lyman-a 
forest gives B H -i <^ 10~ 9 G which is stronger than the limit for a homogeneous universe, 
B h x ^ 2 x 10~ 8 G; while for 50 Mpc coherence length, the inhomogeneous case gives 

-B50Mpc ^ 6 x 10~ 9 G while the homogeneous limit is -B^Mpc & 10 G and for coherence length 
equal to the Jeans length, B\ d <; 10~ 8 G for the Lyman-a case while B% <; 10~ 6 G. 



2. The Electron Density Distribution 

For distant extragalactic sources, the integral in Eq. (1) can be written as an integral over 
redshift using the radial displacement-redshift relation: 

dl( - z) ~ ~ h (i + z)y/n M (i + zf + n A + (1 - n M - n A )(i + zf (2) 

where dl is the proper displacement at redshift z, Hq is the present Hubble parameter, is 
the ratio of the matter density to the critical density in the Universe, and f^A is the vacuum 
energy density over the critical density (see, e.g., Peebles 1993). Following the usual choices for 
the cosmological parameters, we assume h = 0.65, Qm = 1, and = below. (We have also 
considered the case of = 0.2 and ^a = 0.8, which changes the result only by a factor <; 2.) 

The electron density along the line-of-sight to a QSO depends on the mass distribution, the 
baryon fraction, and the ionization fraction. Observations of the Lyman-a forest in QSO spectra 
give evidence of the neutral hydrogen distribution along the line-of-sight which we assume to trace 
the bulk of the baryon distribution. Given a value of the baryon number density over the critical, 
Jib, we can model the electron distribution by assuming an ionization fraction. We assume that 
the bulk of the baryons is ionized hydrogen and helium (24%) with an average density given by 
fib and that the spatial dependence of the baryon distribution (hence, the electron distribution) is 
the same as the neutral hydrogen distribution. 

In modelling the intergalactic medium, we use the analytical approximation of a log-normal 
density distribution for the Lyman-a forest (Bi & Davidsen 1997, and Coles & Jones 1991) which 
gives a simple and consistent description of the observations. The log-normal distribution sets the 
baryon density fluctuations at the Jeans length for photoionized clouds, Aj, which can be written 
in comoving coordinates as (Gnedin & Hui 1996, Bi & Davidsen 1997, and Bond, Szalay, & Silk, 
1988): 

27r / 1 4 K 

~ ~ ■ l ^{l + z)—^hM VC -' . (3) 
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The log-normal distribution gives the electron density distribution as a function of one 
parameter, a(z) 2 , which is the linear variance over the Jeans length. For standard CDM, a good 
fit to a(z) can be written as 

a(z) = 0.08096 + 5.3869(1 + z)' 1 - 4.21123(1 + z)~ 2 + 1.4433(1 + z)~ 3 , (4) 

while ACDM gives a(z) within 20% of standard CDM (Bi & Davidsen, 1997). 

We simulate 10,000 random lines-of-sight to high-redshift QSOs by randomly drawing electron 
densities for cells given by the Jeans length. For each Jeans length cell, we calculate the electron 
over-density (8) with the log-normal equation: 

V ' s/2^a(l + 6) \ 2d 2 J v ' 

The physical electron density is then given in terms of the over-density calculated above, 

n e (z) = n e (z)(l + 8) , (6) 

where the mean electron density as a function of redshift, n e (z) = n e (0)(l + z) 5 , and 
n e (0) = 9.1 x 10~ 6 cm~ 3 n b h 2 . We set fl b h 2 = 0.02 from D/H measurements in high-redshift QSO 
systems (Buries & Tytler 1998a, 1998b) together with the predictions of big bang nucleosynthesis 
(see, e.g., Copi, Schramm, & Turner 1995). This gives n e (0) = 1.8 x 10~ 7 cm~ 3 independent of Hq. 

We add the rotation measures contributed by each Jeans length along the line-of-sight from 
the QSO redshift to z=0. In each cell, the rotation measure, 

where Al = Aj/(1 + z) and can be related to Az through Eq. (2). 

We assume that the large scale magnetic field we are constraining predates the formation of 
structure and is considered frozen into the plasma. The conductivity of the plasma in the Universe 
is large enough for the frozen in assumption to hold on scales ^> 1 A.U. (Cheng & Olinto, 1995) 
which is well below the scales of interest here. The field direction is then taken to be fixed between 
different magnetic field coherence lengths, L c , while the field magnitude changes as a function of 
the electron density. For simplicity, we assume that the overdensities are spherical on average and 
thus scale the field magnitude as B oc n'J^ . 

Therefore, the magnetic field projected along the line-of-sight, B^(z), is fixed within each cell 
to be Si I (2) = B(z) cosO = Bo(n e /n e (0)) 2 / 3 cos9, where Bq is the magnitude of the cosmological 
field in comoving coordinates (actual value at z = 0) and is the angle between the line-of-sight 
and the direction of the cosmological field. is chosen randomly for each magnetic field coherence 
length scale. For magnetic field coherence length, L c = Aj/(1 + z), the random angle of the field 
is chosen at each step, while for L c = N\j, the random direction is changed after N steps. 
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The coherence length for a given cosmological field depends on the particular magnetogenesis 
model and the evolution of the field up to recombination. Given the damping of the short 
wavelength modes up to the Silk scale (Jedamzik, Katalinic, & Olinto 1998) the coherence scale 
is likely to be above A j, which today is A j ~ 1 Mpc. We chose to study three choices of L c that 
span the range of interest for plausible cosmological fields: Aj, 50 Mpc, and Hq 1 . As discussed 
below, our results are not very sensitive to the coherence length of the field and our conclusions 
can be applied to any model of magnetogenesis where L c ^ Aj. 

3. Results 

After integrating Eq. (1) for 10,000 lines-of-sight per 0.1 redshift interval, we display our 
results in Figs. 1 and 2, where we plot the median of the absolute value of the RM {median |RM|, 
circles in Fig. 1), the mean absolute value of RM (< |RM| >, squares in Fig. 1), and the variance 
(< RM 2 >, triangles in Fig. 2). For comparison, we also show the same functions for the case of 
a homogeneous electron distribution throughout the universe {median |RM| are the solid lines in 
Fig. 1, the < |RM| > are the dashed lines in Fig. 1, and the variance < RM 2 > is the solid line in 
Fig. 2). We also bin the available data on QSO rotation measures in two redshift bins, z = — 1.5 
and z = 1.5 — 2.5 and display the < |RM| > (dashed-dotted line) and the median\RM\ (dotted 
line) in Fig. 1. 

Not surprisingly we find that the Lyman-a forest introduces a very large scatter in the RM 
for different lines-of-sight. Not only the magnitude of the variance is many orders of magnitude 
greater than the case of a homogeneous universe but the variance and the mean RM vary 
significantly within our unrealistically large simulated sample of 10,000 lines-of-sight per Az = 0.1 
redshift interval (Fig. 2). With such a large variance, the limited data set available at present 
is not sufficient to constrain well the cosmological magnetic field. In fact, Oren & Wolfe (1995) 
showed how sensitive is the redshift dependence of the variance to different binning procedures 
and how the observed distribution is non-gaussian. The trends in the variance as a function of 
redshift in Fig. 10 of Oren & Wolfe (1995) (where they plot different binning of data compiled by 
Welter, Perry, & Kronberg 1984) are incompatible with the variance behaviour in a homogeneous 
universe (see Fig. 2). In particular, the homogeneous case has a monotonically increasing variance 
with redshift while the inhomogeneous case can reproduce the behaviour shown in their Fig. 10. 
In this case, the non-gaussian tail of the RM is due to a few large density clouds along random 
lines-of-sight. The non-gaussian tail can also be due to non cosmological magnetic fields such as 
small scale structure unaccounted for in the Galactic magnetic field. 

Of the statistical indicators we studied, the variance is the noisiest followed by the mean and 
the median is the most stable indicator. This behaviour is evident in Figs. 1 and 2. The fact that 
the mean and the mediam differ significantly in the case of an inhomogeneous universe is further 
evidence of the non-gaussian behaviour of the rotation measures. The mild magnetic coherence 
length dependence of the median and the mean, is such that as the coherence length grows the two 
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indicators approach each other. We use the best statistical indicator, the median\RM\, to set our 
limits for both the homogeneous and inhomogeneous cases. We only set limits versus detections 
due to the presently unsurmountable challenge of separating locally generated contributions to the 
rotation measure from cosmologically generated ones. 

The available date on QSO RMs is sparse and dominated by the Galactic magnetic field 
contribution. Courageous attempts at modelling and subtracting the Galactic contribution can 
be found in Simard-Normandin & Kronberg 1980, Simard-Normandin, Kronberg, k, Burton 1981, 
Kronberg & Perry 19982, Welter, Perry, & Kronberg 1984, Watson & Perry 1991, Vallee 1990, 
1993, and Oren & Wolfe 1995. (For a critical comparison of these attempts see Oren & Wolfe 
1995.) Given the sparseness of the available data, we chose to bin the Welter, Perry, Sz Kronberg 
1984 catalog in two redshift bins (0 < z < 1.5 and 1.5 < z < 2.5) and calculated the mean absolute 
RM (dashed-dotted line) and the median absolute RM (dotted line) for this data and compared 
these with our predictions. 

For the £Im = 1 hat universe (with n b h 2 = 0.02) we get for cosmological fields coherent 
accross the present horizon, B H -i <^ 10~ 9 G in the case of a Lyman-a forest which is stronger than 
the limit for a homogeneous universe, B h 1 ^2x 10~ 8 G; while for 50 Mpc coherence length, the 

inhomogeneous case gives i?50Mpc ^ 6 x 10 -9 G while the homogeneous limit is B^ 0Mpc <; 10~ 7 G 
and for coherence length equal to the Jeans length, B\j <J 10~ 8 G for the Lyman-a case while 
B\j & 10~ 6 G. Changing cosmological parameters to, for example, = 0.2, $7 a = 0.8 in a flat 
universe, the results only change by a factor of 1.6; we find median \RM\ (Qm = 0.2,$1\ = 0.8) 
~ 1.6 x median\RM\ {Qm = 1>^A = 0), and the limits on cosmological magnetic fields is 1.6 
times lower. These limits are slightly stronger than limits attained by alternative methods. For 
instance, the limit due to the isotropy of the cosmic background radiation is B R -i <; 4.410~ 9 G 
for VLq = \ and h = 0.65 (Barrow, Ferreira, Silk 1997). Given the non-gaussian nature of the RM 
distribution, the sparseness of the data, and the difficult Galactic subtraction, it is encouraging to 
find the RM limit for the largest coherence length to be consistent with and comparable to the 
cosmic background radiation limit. 

Our results imply stronger constraints on scenarios of primordial magnetic field generation 
and evolution by limiting the spectrum of primordial fields today both the magnitude and the 
scale dependence, i.e., we limit Bq(L c ). This lowers the available primordial seed for galactic 
dynamos and for structure formation by primordial fields. Our results also help the prospective 
of detectors pointing back to ultra-high energy cosmic ray (UHECR) sources, by limiting the 
primordial component of the intergalactic medium field. Magnetic effects on UHECRs will be 
more important if post recombination fields, limited to sites of structure formation, are much 
stronger than the fossil fields from the early universe (e.g., Ryu, Kang, & Biermann 1998). 
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a fellowship of the Istituto Nazionale di Fisica Nucleare at the University of Chicago. A.V.O. was 
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Fig. 1. — Observed and predicted statistics of QSO rotation measures for magnetic field coherence 
lengths of the Jeans length, 50 Mpc, and Hg 1 . The predicted absolute mean (squares) and median 
(circles) rotation measure for a 6 nG field in the Ly-a forest approximation drawn from 10000 lines 
of sight at each redshift. The continuous lines represent the predicted absolute mean (solid line) and 
median (dashed) for a homogeneous electron density and a 6 nG field. The two bins (0 < z < 1.5 
and 1.5 < z < 2.5) show the observed absolute mean RM (dashed-dotted) and median (dotted) in 
the dataset of WPK. 

Fig. 2. — The variance of QSO rotation measures predicted for the Ly-a forest (triangles) and a 
homogenenous density (solid line) for a 6 nG field and the coherence lengths shown in Figure 1. 
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